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Abstract. In an incompressible velocity field, the surface area of a vol- 
ume varies with time, but volume remains unchanged. If incidentally the 
surface becomes spherical along time, the area reaches a local minimum, 
since sphere has the least area that surrounds a volume. So the area is 
a function of time that is locally convex at this point. When applied to 
an incompressible Navier-Stokes fluid, this property is used to compute 
an inequality that suggest a criterion to non-existence of initial configu- 
rations of velocity fields, revealing its impossibility to evolve with time. 
Three velocity fields are proposed as examples. One of them agrees the 
inequality, the other two violates it. 

Keywords. Spherical Surface Area, Isoperimetric Inequality, Navier-Stokes 
Equations. 



1. Introduction 

Every dynamical system, described by differential equations, deals with the 
initial value problem. This is, given an initial condition, one tries to deter- 
mine whether the system can evolve with time or not beginning from that 
condition. Sometimes, it may be possible to determine whether there are one 
or several solutions for equations with the initial condition. There was pro- 
posed in jT] that initial conditions for incompressible Navier-Stokes velocity 
fields are useful to find its time evolution, in such a way that given suitable 
restrictions to the initial velocity field, the system is determined at least for 
any finite time after. In the same way, Beale, Kato and Majda [2] proved that 
a smooth velocity field may lose its regularity some time after, in such a way 
that the maximum vorticity becomes unbounded. Hence, to find properties 
of the initial velocity field is a challenge. In this paper we propose a criterion 
to the non-existence of some of these velocity fields. 
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2. Transport theorem for surfaces 

Reylods transport theorem |3j is a very useful tool since it lets introduce the 
time derivative of a volume integral inside the integrand of a static volume 
integral. We would like to do the same with a surface integral. It is, to trans- 
form the time derivative of a surface integral, which surface is moving and 
changing its shape, and to obtain a fixed surface integral with a time deriva- 
tive inside its integrand. The next theorem shows how to find this issue (see 

a)- 

Theorem 2.1. Let u(t,x) £ R 3 be a velocity field with components Ui that are 
enough smooth, and let be f (t, x) € R also a smooth function. Let fl C R 
be a region of the field with boundary dCl. The unitary normal vector to dfl 
is ft, with components n,. Then, 



d_ 
dt 



fd 2 x = / [d t f + mdj + (dm - eijmrij) f] d 2 x, (2.1) 



where eij = \ (diUj + djUi) is the infinitesimal strain tensor (defined by, e. 
g., W- 

Proof. Since the vector normal to the surface is unitary and the surface 
closed, we can use Gauss theorem 

— / fd 2 x — — ( nifnid 2 x 
dt Jdn{t) dt y 9 Q( t ) 



d 
dt 



8, (fm) d 3 x. (2.2) 

Q(t) 



Then, we can apply Reynolds transport theorem to the volume integral 
d 



dt 



fdx=j [dtdi{fn i )+u j d,jd l {fn i ) + djUjd l {fn i )]d i x. (2.3) 

an(t) 



Using the chain rule twice in the second term in the integral of right hand 
side and taking in to account that time and space derivatives commutes in 
the first term of right hand side, 



d 
dt 



fd x= I {dt [(d t + ujdj) (frii)} - dj (fn^Uj)} d x + 



+ / d i (fn i d jUj )d 6 x. (2.4) 



So, we can use again Gauss theorem, 



/ fd 2 x = I f (d t + Ujdj) \ \tiini ) d 2 x 
an(t) Jon \^ 



dt 



+ / {niUi [(d t + Ujd 3 ) f + djUjf] - fdiUjUjUi} d x. (2.5) 
Jan 
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The first term of right hand side is the derivative of a constant and it vanishes. 
Then, from 

d_ 
d~t 



fd 2 x 



d€l{t) 



Oil 



\(d t + Uidi) f + (diUi - d t u j n J n i ) f] d 2 x, 



the relation (2.1 1 arises and the theorem is proved. 



(2.6) 
□ 



Equation (2.1| is similar to the transport theorem for moving surfaces 
[6|- [7], which usually is written in terms of both, normal velocity and cur- 
vature of the surface. Now that we know the rate of change of the surface 
integral of a magnitude with time, we would like to know whether the area of 
the surface grows, diminishes or remains constant with time when the volume 
does not change. We knows a particular case yet. One of the properties of 
the sphere is that it has the least area that encloses a volume. So, the area 
of the sphere only can increase or be the same few time after. This means 
that the area is a convex function of time near of the minimum. The next 
theorem refletcs this situation. 

Theorem 2.2. Let u(t,x) £ R 3 be a velocity field with components u r ,ug,u,p 
in spherical coordinates. Let S 3 C R 3 be a spherical region of the field with 
boundary S 2 and radium r. Also, there exists only one region fl (t) C R 3 for 
t ^ to such as Q(t) — > S 3 when t — » to. For each S 3 , and every t, if the 
velocity field holds the incompressibility statement, V • it = 0, then 



p7T p2lT 


' 2 D 


'o Jo 





r 2 sin 



> 0, 



(2.7) 



where e r 



d r u r ( for stain tensor in spherical coordinates see [HJj. 



Proof. Taking into account the very well known isoperimetric inequality for 
three dimensions [9]-[10], we have 



d 2 x > 3 



d A x 



n(t) 



(2.8) 



where the equality holds for the sphere S 3 . We substract the area of S 2 on 
both sides, 



d 2 x 



d 2 x > 3 



>3, r 









f 


f d 3 x 


- [ d 2 x 




Jn(t) 


Js 2 



d 3 3 



/n(t) 



(2.9) 



Due to the incompressibility of the fluid, S 3 and ft have the same volume. 



The right hand side of (2.9) then vanishes 



d 2 . 



d 2 x > 0. 



s 2 



(2.10) 
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In addition, the area time derivative is given by (2.1 ), with / = 1 and diUi — 
0, 



d r u r r sin 



' d [ / 




— ax 
at Jdn(t) 


=- r i 


Jo Jo 



-d r 



u r r sin 



-d r 



diUid x 



JO 

2 



diUid x = 0. 



S 3 



(2.11) 



So the area of a sphere reaches its minimum at time t — t in a incompressible 



velocity field. This property together with (2.10) means that the area is a 
local convex function of time in a range close to to- Therefore, the second 
time derivative of this function at to holds 



dt 2 



d 2 :, 



ao(t) 



(to) > 0. 



(2.12) 



The second time derivative of the area can be computed applying (2.1 ) twice 

(—riinjdiUj) d 2 x 



d?_ 
dt 2 



d 2 :, 



ao(t) 



(to) 



(niUjdiUj 



dQ(t) 

D 
Dt 



{riinjdiUj) > d 2 x 



D 



(to) 
(to) 



(2.13) 

where in the last line we have used that the normal vector to the surface of 



(d r u r ) - — (d r u r ) } d x, 

S 2 I Lit 



the sphere only has radial component. Taken (2.12) together with (2.13), we 



can find (2.7) at time t = to. But the spherical surface is independent of time, 



since the time dependency is in the integrand of the last line of (2.13). This 



is, at every time, for every spherical surface, there exist a volume, which is 

;ry 
□ 



a function of time, that converges to the sphere. Then (2.7) is held at every 
instant of time. 



Given that we have a surface integral, it does not matter what is the 
velocity distribution inside the sphere but just that velocity distribution on 
its surface. Therefore, this theorem asserts that if there exist at least a sphere 



in the domain of the incompressible velocity field that violates ( 2.7 ), evolution 
with time is forbidden for that velocity field. The next lemma applies this 
theorem to incompressible Navier-Stokes fluids. 

Lemma 2.3. Let u(t,x) G R 3 be an incompressible velocity field, V • u = 0, 
which evolves in time according to the Navier-Stokes equations 

d t u + u-Vu = ^Au-Vp. (2.14) 
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Here, p is the pressure, the density is p = 1 and /i is the viscosity. Veloc- 
ity components in spherical coordinates are denoted by u r ,ug,u,p (9 and <j> 
are polar and azimuth angles, respectively). Then, at every time t, for every 
spherical region of the field S 3 C R, 3 with boundary S 2 and radius r, we have 



s 2 



jd 2 p(i, x) - F{t, x) - fj,G(t, x) + {d r u r (t, f)) 2 } d 2 x > (2.15) 



where 



F(t,r,d,< 



(d r u r ) 2 — d r ) dgu r 



-d r (?*) ^. 

\ r J si 



sin 9 



G(t, r, 6, <j>) = dlu r + 2d 2 (— ) + (dp r Ur + cot 9d e d r u r ) 



(dgU r + COt 9dgU r ) 



(2.16) 



(2.17) 



Proof. The radial direction of the equation (2.14) is given by 



d t u r + u r d r u r H <9git r 

r 



-<9 r p + p 



<9 2 u r 



2<9 r 



r sin# 
1 



u 



^d 2 9 U r 



cot I 



d e u r 



(2.18) 



(see, e.g., [TT])Next we can take the partial derivative of this relation with 
respect to r, and then we use the identity 



-prr = o t Q + u r d r Q H <9 e Q H ^- 

Vt r r sm t 



(where Q is a scalar magnitude) to group terms, obtaining 



D 
Dt 



(d r u r ) = -d 2 P + F + nG, 



Substitution of this relation on (2.7) gives rise to (2.15). 



(2.19) 

(2.20) 
□ 



This lemma means that if we find at least a spherical surface for which 
the incompressible velocity field does not hold (2.15 ), that field can not evolve 
acording to Navier-Stokes equations. Notice that the lemma is only useful 
when the inequality is violated. Lets see it with three examples. 

Example (1). At time to, let a velocity field be given by 



(2.21) 



(in spherical coordinates) inside a bigger sphere of ratio R. The fluid of this 
velocity field spins around the z axis and is divergent-free. We would like to 
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Figure 1. Coordinates origin shifted a distance L along y axis. 



confirm that (2.15) is correct. The computation of (2.16) and (2.17) to this 



velocity field gives 



F{r,6, 
G{r,6, 



sin 2 (9 




(2.22) 
(2.23) 



Computation of double radium derivative of pressure is more difficult. We 
can work out the pressure, as usual, by solving the Poisson equation ob- 
tained when we take the divergence of incompressible Navier-Stokes equations 



(2.14). So, this non-local function of spatial derivatives of velocity is 

V ■ ( u ■ Vu) 



V (t, #) = 



1 

4tt 



d 3 3 



R3 



(2.24) 



In our case, the corresponding derivations and integration in the sphere of 
radium R gives us 

p(r,M)= \ r2 - R2 ' ( 2 - 25 ) 

being R > r. Then, the double time derivative of the surface area that con- 
verges to a sphere of radium r in this velocity field is 

J \d 2 rP - F - fiG+ {d r u r f}r 2 sin 9d9d<f> 



-- sin 2 9\r 2 sm6d6d<j) = 0. 



(2.26) 



This time, the result does agree with the inequality (2.15). But still remains 
regions of the R 3 where we can look for violation of the inequality. Now, the 
system of reference is shifted a distance L from the z axis along y axis (see 
Fig. [l]), instead of be on it. With the identities given by 



r cos 9 = r' cos 9' 

r sin 9 cos <j> = r' sin 9' cos <j>' 

L = r' sin 9' sin <b' — r sin 9 sin < 



(2.27) 



we can change 
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Up 



l_/2 
3' 



R 2 



(2.28) 



by 



u r — L sin 9 cos 
ue = L cos cos 
M0 = r sin 



L sin ( 



- 2 + L 2 + 2Lrsin( 



R 2 



(2.29) 



with R > L + r. Of course, this velocity field still has V • u = 0. Repeating 
again the steps like before, we find 



Ffr, 



G(r,9,< 



L 



L 2 



L cos < 

t3 Qin 



! sin 2 9 — cos 2 ( 



(sin 2 4> + cos cos 4> sin 0()2.3O) 
(2.31) 



and hence, 



71 />27T 



Jo 



{d 2 p-F-fiG+(d r u r f} 



r 2 sin 



(2.32) 



This result also agrees with the inequality (2.15|. Moreover, it is independent 
of the radium R. So, this inequality can be extrapolated to R 3 doing R — >• 00 
or L — > 00 and using revolution symmetry around z' axis. 



Example (2). In this example, we will see that the divergent-free velocity field 
(see Fig. [2]) given in spherical coordinates, at a time to ,by 



u r = 

u g = (2.33) 
i*0 = r k sin#, 



with 2 < fc 6 N, does not hold the inequality (2.15) for every sphere that is 
inside a bigger sphere of radium R. As before, first we compute F, G from 
velocity and its derivatives, and then, p from them and from the integral over 
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the sphere of radium R (being R> r). So, we obtain 



F(r,0,4>) 
G(r,9,cf>) 

p{t,6,4>) 



(2k- l)r 




2fc - 2 sin 2 « 



(2.34) 
(2.35) 



E 

nGN°\{2fc} 



R 



2fc-r 



2n + 1 



r 2k-n 2k 



„2k 



1 J 2k 



f [(k - 1) sin 2 6' + 1] sin 0'P n (cos (0-0')) dO' 
Jo 



— ^ In — 

2A: + 1 r 



,2 A' 



/ [(Jfe - 1) sin 2 0' + 1] sin 0'P 2k (cos (0 - 0')) dO' 
Jo 



(2.36) 

(where N°\ {2k} = 0, 1, 2. ..2k— 1, 2fc+l, ... and P n (x) are Legrendre polimom- 
inals) and hence, 



7T /•2-7T 



o Jo 



{<9 2 p - F - fj,G+ (d r u r ) 2 } r 2 sinOd0d<j) 



-2irr 2k 

n6N°\{2fc} 

-2irr 2k 



2k(2k-l) 2k (2k -I) n(n-l)R 



2k-nl 



2k + n+l 
2k (4k - 1) - 1 



2k — n 2k — n r 2k ™ 

R 



4fc + 1 



2k (2k- l)ln 



A 2k - (2k - 1) -Trr 



2/,' 



(2.37) 



where 



/•7T /*7T 

A n = / (sin 2 6'+l)sin6lsin6l'P„(cos(6l-6''))d^'d6 , ' (2.38) 
Jo Jo 

Notice that, since —2 < A n < 2 and r < R, the summatory converges. 
Moreover, we have used R as a parameter to compute the pressure and we 
can make it very large. When R ^> r, we can approach (2.371 by 



-i 7->2fc-2 
^ 2k . ■ afc _2 ^ + Q 



fc — 1 r 2fc 



(2.39) 



However, since = 12/5, it is impossible that it holds the inequality 



2irr 



2k 



1 R 



2k-2 



k — 1 r 2k 2 



a 2 + o 



7? 



2fc-3 



> 0. 



(2.40) 



So we conclude that surprisingly the velocity field (2.33) can not evolve ac- 
cording to incompressible Navier-Stokes equations. 
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Figure 2. Velocity profile of (2.33) in the x-y plane 



Example (3). In this last example, we will see that the divergent-free velocity 
field (see Fig. [3]) given, at a time to ,by 



r sin v cos 



u r = (R 2 - 

ug = — (R 2 — r 2 sin^ 9) s'mt 
= 0, 



(2.41) 



with < r < oo and < 9 < tt, does not hold the inequality (2.15) for every 



sphere. Proceeding as before, we again compute F, G and p, this one worked 
out throught an integral over all three dimensional space. So, we obtain 



F (r, 6, 4>) = (r 2 sin 2 9 - R 2 ) (r 2 sin 2 9 - cos 2 9) sin 2 
G(r,0,<P) - ( 1 °S-2sin 2 9^ COS9 



p(r,6,<t>) = 0, 



(2.42) 
(2.43) 

(2.44) 



and hence, 

«7T p27T 



\d 2 rP -F-iiG+ (<9 r w,.) 2 }r 2 sir 



I6nr 2 



R 2 



(2.45) 



However, the inequality (2.15) does not hold when the radius of the probe 
sphere is 

77- 



r>~R. 



(2.46) 



Then, the velocity field (2.41) can not evolve according to incompressible 
Navier-Stokes equations. 
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FIGURE 3. Velocity profile of (2.41 ) in a plane that contains 
the z axis. 



3. Conclusion 

We have shown that given an incompressible velocity field at a initial time, we 
can test whether its time evolution is forbidden by a criterion. It is related 
with the non-negativeness of the double time derivative of the area of a 
volume that becomes a sphere at that instant. Of course, if velocity agrees 
the inequality at a give time it also agrees that the reminder of the time, 
because in other case, the field had not evolved to reach that instant. In 
particular we have worked the inequality out to a Navier-Stokes fluids. We 
also have found two particular incompressible velocity fields that can not 
evolve acording to Navier-Stokes equations. 
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